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1 Introduction 



The charmless B decays arouse more and more interests recently, since it is a good place 
for study of CP violation and it is also sensitive to new physics [0. Factorization approach 
(FA) is applied to hadronic B decays and is generalized to decay modes that are classified 
in the spin of final states p|, ^, ^. FA gives predictions in terms of form factors and decay 
constants. Although the predictions of branching ratios agree well with experiments in most 
cases, there are still some theoretical points unclear. First, it relies strongly on the form 
factors, which cannot be calculated by FA itself. Secondly, the generalized FA shows that the 
non-factorizable contributions are important in a group of channels p, The reason of this 
large non-factorizable contribution needs more theoretical studies. Thirdly, the strong phase, 
which is important for the CP violation prediction, is quite sensitive to the internal gluon 
momentum . This gluon momentum is the sum of momenta of two quarks, which go into two 
different mesons. It is difficult to define exactly in the FA approach. To improve the theoretical 
predictions of the non-leptonic B decays, we try to improve the factorization approach, and 
explain the size of the non-factorizable contributions in a new approach. 

We shall take a specific channel B ^ mi as an example. The B ^ mi decays are responsible 
for the determination of the angle 02 in the unitarity triangle which have been studied in the 
factorization approach in detail 0, ^, The recent measurements of -B ^ tt+tt^ by CLEO 
Collaboration attracted much attention for these kind of decays . The most recent theoretical 
study [0 attempted to compute the non-factorizable diagrams directly. But it could not also 
predict the transition form factors of B ^ n. 

In this paper, we would like to study the B ^ mr decays in the perturbative QCD approach 
(PQCD) [^. In the B ^ mi decays, the B meson is heavy, sitting at rest. It decays into two 
light mesons with large momenta. Therefore the light mesons are moving very fast in the rest 
frame of B meson. In this case, the short distance hard process dominates the decay amplitude. 
We shall demonstrate that the soft final state interaction is not important, since there is not 
enough time for the pions to exchange soft gluons. This makes the perturbative QCD approach 
applicable. With the final pions moving very fast, there must be a hard gluon to kick the light 
spectator quark d or u (almost at rest) in the B meson to form a fast moving pion (see Figure |lD. 
So the dominant diagram in this theoretical picture is that one hard gluon from the spectator 
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quark connecting with the other quarks in the four quark operator of the weak interaction. 
Unhke the usual FA, where the spectator quark does not participate in the decay process in a 
major way, the hard part of the PQCD approach consists of six quarks rather than four. We 
thus call it six-quark operators or six-quark effective theory. Applying the six-quark effective 
theory to B meson exclusive decays, we need meson wave functions for the hadronization of 
quarks into mesons. Separating that nonperturbative dynamics from the hard one, the decay 
amplitudes can be calculated in PQCD easily. Most of the nonperturbative dynamics are 
included in the meson wave functions, but in the correction that soft gluon straddle the six- 
quark operators, there are some nonfactorizable soft gluon effects not to be absorbed into the 
meson wave functions. Such effects can be safely neglected in the B meson decays [Q. 

Li performed the calculation of B^ tt+vt" in ref. [|r^ using the PQCD formalism, where the 



factorizable tree diagrams were calculated and the branching ratios were predicted. In another 



paper [^], Dahm, Jakob and Kroll performed a more complete calculation, including the non- 
factorizable annihilation topology and the three decay channels of B ^ tttt decays. However, the 
predicted branching ratios are about one order smaller than the current experiments by CLEO 
B. In connection with this, Feldmann and Kroll concluded that perturbative contributions to 



the B ^ n transition form factor were much smaller than nonperturbative ones [12|. As we 



shall show later, the pion wave function must be consistent with chiral symmetry relation 

- q^{0\u-ff,-f5d{x)\7r'{q)) = {m.^ + md)(0|M75rf(x)|7r"(g)). (1) 

This introduces terms that were not considered in above calculations. In this paper, considering 
the terms needed from chiral symmetry, we calculate the B ^ n transition form factors and 
also the non-factorizable contributions in PQCD approach. We then show that our result for 
the branching ratio B tt+tt^ agree with the measurement. Among the new terms, it is 
worthwhile emphasizing the presence of annihilation diagrams which are ignored in FA. We 
find that these diagrams can not be ignored, and furthermore they contribute to large final 
state interaction phase. 
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2 The Frame Work 



The three scale PQCD factorization theorem has been developed for non-leptonic heavy meson 
decays |jl3[, based on the formalism by Brodsky and Lepage |]14|, and Botts and Sterman |jl5 



The QCD corrections to the four quark operators are usually summed by the renormalization 



group equation |jT6[. This has already been done to the leading logarithm and next-to-leading 
order for years. Since the b quark decay scale rub is much smaller than the electroweak scale mw, 
the QCD corrections are non-negligible. The third scale 1/6 involved in the B meson exclusive 
decays is usually called the factorization scale, with b the conjugate variable of parton transverse 
momenta. The dynamics below 1/b scale is regarded as being completely nonperturbative, and 
can be parametrized into meson wave functions. The meson wave functions are not calculable 
in PQCD. But they are universal, channel independent. We can determine it from experiments, 
and it is constrained by QCD sum rules and Lattice QCD calculations. Above the scale 1/b, the 
physics is channel dependent. We can use perturbation theory to calculate channel by channel. 

Besides the hard gluon exchange with the spectator quark, the soft gluon exchanges between 
quark lines give out the double logarithms ln^(P6) from the overlap of collinear and soft diver- 
gences, P being the dominant light-cone component of a meson momentum. The resummation 
of these double logarithms leads to a Sudakov form factor exp[— s(P, 6)], which suppresses the 
long distance contributions in the large b region, and vanishes as 6 > 1/Aqcd- This form factor 
is given to sum the leading order soft gluon exchanges between the hard part and the wave 
functions of mesons. So this term includes the double infrared logarithms. The expression of 
s{Q, b) is concretely given in appendix |B[ Figure ^ shows that e~'^ falls off quickly in the large b, 
or long-distance, region, giving so-called Sudakov suppression. This Sudakov factor practically 
makes PQCD approach applicable. For the detailed derivation of the Sudakov form factors, see 
ref-BH. 

With all the large logarithms resummed, the remaining finite contributions are absorbed 
into a perturbative b quark decay subamplitude H{t). Therefore the three scale factorization 
formula is given by the typical expression, 

'■^ dp. 

'1/b fl 

where C{t) are the corresponding Wilson coefficients, <l>(x) are the meson wave functions and 



C{t) X H{t) X X exp 



Jl/b u 



(2) 
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the variable t denotes the largest mass scale of hard process if, that is, six-quark effective 
theory. The quark anomalous dimension 7^ = — cts/vr describes the evolution from scale t to 
1/6. Since logarithm corrections have been summed by renormalization group equations, the 
above factorization formula does not depend on the renormalization scale explicitly. 

The three scale factorization theorem in eq.(0) is discussed by Li et al. in detail Below 
section ^, we shall give the factorization formulae for B ^ tttt decay amplitudes by calculating 
the hard part H{t), channel dependent in PQCD. We shall also approximate H there by the 
0{as) expression, which makes sense if perturbative contributions indeed dominate. 

In the resummation procedures, the B meson is treated as a heavy-light system. The wave 
function is defined as 

*B = + ^Bh5^B{ki, h), (3) 

where A^^c = 3 is color's degree of freedom and (psiki, kb) is the distribution function of the 
4- momenta of the light quark (ki) and b quark {kh) 

Mkuh) = ^^-l=|^e^'=-^(0|T[%)^756(0)]|i?(Pi.)). (4) 

Note that we use the same distribution function 0_b(/ci, kb) for the term and the ms term 
from heavy quark effective theory. For the hard part calculations in the next section, we use 
the approximation nif, c:^ itlb, which is the same order approximation neglecting higher twist 
of {rriB — nibj/mB- To form a bound state of B meson, the condition k^ = Pb — ki is required. 
So (pB is actually a function of ki only. Through out this paper, we take p^ = (p° ±p3)/72, 
Pt — {P^iP^) the light-cone coordinates to write the four momentum. We consider the B 
meson at rest, then that momentum is Pb = (mB/v^)(l, 1,0^). The momentum of the light 
valence quark is written as (/cj*", fcf , ki^), where the ki^- is a small transverse momentum. It is 
difficult to define the function 05(^5*", fcf, ki^-). However, in the next section, we will see that 
the hard part is always independent of k^ , if we make some approximations. This means that 
ki can be integrated out in eq.(|D, the function (pBikf, A;f , ki^) can be simplified to 



()B(xi,kir) = Pb J dkt(j)B{kt,k^ ,kiT) 



pl2V2Nj (27r)3 

where Xi = k^ /p^ is the momentum fraction. Therefore, in the perturbative calculations, we 



do not need the information of all four momentum ki. The above integration can be done only 
when the hard part of the subprocess is independent of the variable kf. 

The 71 meson is treated as a light-light system. At the B meson rest frame, pion is moving 
very fast. We define the momentum of the pion which contain the spectator light quark as 
P2 = (mB/v/2)(l,0,OT). The other pion which moving to the inverse direction, then has 
momentum P3 = (mB/v^)(0, 1, Ot). The light spectator quark moving with the pion (with 
momentum P2), has a momentum (/c^, 0, k2T)- The momentum of the other valence quark in 
this pion is then {P2' — k2, 0, — k2T)- If we define the momentum fraction as X2 = k2 / P2 1 then 
the wave function of pion can be written as 

= 7^^^ [^7r0^(a;2, kar) + mo0U^2, ksr)] , (6) 
where (pn(yX2, k2T) is defined in analogue to eq.(|, ^ and 0^(x2, k2T) is defined by 

0:(x2,k2T) = ^^/^^^^e^(-^^"^--'^--^-)(0|T[rf(0)7^^ (7) 



Note that as you shall see below, ttlq given as 
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mo = (8) 

in eq. (H) is not the pion mass. Since this mo is estimated around 1 ~ 2 GeV using the quark 
masses predicted from lattice simulations, one may guess contributions of mo term cannot be 
neglected because of m-o^ m^. In fact, we will show this mo plays important roles to predict 
the 5 — s> vrvr branching ratios in section |. 

The normalization of wave functions is determined by meson's decay constant 

(0|J(0)7^75w(0)|vr(p)) = «p^^. (9) 

Using this relation, the normalization of 0,^ is defined as 

J rfx2rf^k2r07r(a;2,k2r) = 2^/^ " (-^^^ 
Moreover, from eq.(P) you can readily derive 



m? 



(0M(0)75w(0)|7r(j9)) = -t (11) 

so defining mo such as eq.(^), the normalization of (f)'^ is the same one to eq.(TU). 



The transverse momentum is usually conveniently converted to the b parameter by 
Fourier transformation. The initial conditions of (f)\'\x), i = B, re, are of nonperturbative 
origin, satisfying the normalization 

0f)(x,6 = O)rfx=^-|=, (12) 
with fi the meson decay constants. 



3 Perturbative Calculations 

With the above brief discussion, the only thing left is to compute H for each diagram. There 
are altogether 8 diagrams contributing to the B ^ nn decays, which are shown in Figure |^. 
They are the lowest order diagrams. In fact the diagrams without hard gluon exchange between 
the spectator quark and other quarks are suppressed by the wave functions. The reason is that 
the light quark in B meson is almost at rest. If there is no large momentum exchange with 
other quarks, it carries almost zero momentum in the fast moving vr, that is the end point 
of pion wave function. In the next section, we will see that the pion wave function at the 
zero point is always zero. The Sudakov form factor suppresses the large number of soft gluons 
exchange to transfer large momentum. It is already shown that the hard gluon is really hard 
in the numerical calculations of B Kn |TB|. The value of a^/vr is peaked below 0.2. And in 
our following calculation of B —>■ decays this is also proved. 

Let's first calculate the usual factorizable diagrams (a) and (b). The four quark operators 
indicated by a cross in the diagrams, are shown in the appendix A. There are two kinds of 
operators. Operators Oi, O2, O3, O4, O9, and Oio are (V — A){V — A) currents, the sum of 
their amplitudes is given as 

rl roc 

Fe = —IQuCpTn^B I dxidx2 / bidbib2db2 4'Bixi,bi) 
Jo Jo 

X {[(1 + X2)(l)n{x2, b2) + (1 - 2x2)0Ua;2, &2)r^] 

xas{tl)K{xi,X2, bi, 62) expl-Ssitl) - Sl{tl)] + 2r^0^(x2, ^2) 

X asitl)Kix2, xi, 62, bi) exp[-SBitl) - SUtl)]] , (13) 

where r^^ = mQ/rriB = m'^/['mB{iTLu + 171^)]. Cp = 4/3 is a color factor. The function 
he{xi,X2,bi,b2) and the Sudakov form factors SB(ti) and S-^iti) are given in the appendix 
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^. The operators O5, Oe, O7, and Og have a structure of (V — A)(y + A). The sum of their 
amphtudes is 

Ff = — 327rCirm|r7r / dxidx2 / bidbib2db2 (f)Bixi,bi) 

Jo Jo 

X {[07r(a;2,&2) + (2 + X2)0^(x2,&2)^7r] (t^) /ig (Xi , , &1 , ^2) 

X exp[-5B(te) - '^^(^e)] + [3^10.(0:2, ^2) + 2(1 - a;i)0;(x2, 62)r.] 

X a,(t2)/i^(a;2, xi, 62, &i) exp[-SB(t2) - Sl{tl)]} . (14) 

They are proportional to the factor r^. There are also factorizable annihilation diagrams (g) 
and (h), where the B meson can be factored out. For the {V — A)(y — A) operators, their 
contributions always cancel between diagram (g) and (h). But for the {V — A){V + A) operators, 
their contributions are sum of diagram (g) and (h). 



— 647rCi?m|r^ / dx2dxz 
Jo 



b2db2b3db3as(ta)ha{x2, xs, 62, ^3) 



X [2<f)^{x2, b2)<j)'^{x3, 63) + X2MX3, h)<P'^{x2, 62)] exp[-Sl{Q - SliQ] , (15) 

These two diagrams can be cut in the middle of the diagrams. They provide the main strong 
phase for non-leptonic B decays. Note that vanishes in the limit of mo = 0. So the mo term 
in the pion wave function does not only have much effect on the branching ratios, but also the 
CP asymmetries. Besides the factorizable diagrams, we can also calculate the non-factorizable 
diagrams (c) and (d) and also the non-factorizable annihilation diagrams (e) and (f). In this 
case, the amplitudes involve all the three meson wave functions. The integration over 63 can be 
performed easily using 6 function 5(63 — 61) in diagram (c,d) and 5(63 — 62) for diagram (e,f). 



32 I /"i r°° 

Me = —-n:CF\2Ncmi dxidx2dx3 / 6i(i6i62C?&2 0B(a;i, 6i)07r(a;2, ^2) 

6 ^ Jo Jo 



32 

3 ' " Jo Jo 

x0^(x3, bi)x2as{td)hd{xi,X2, X3, 61, 62) exp[-S'ij(td) - Sl{td) - Sl{td)] , (16) 

32 I r°° 

M.^ = — TiCFy2Ncm\ / dxidx2dx3 / 6i(i6i62rf&2 ^^(a^i, &i) 
3 * JO Jo 

X |- X2(t)n{x2, b2)(j)n{x3, 62) + {X2 + X3)(f)'^{x2, b2)(f)'^{x3, 62)^^ 

X as{t})hf\xi,X2,xs,bub2)exp [-^^(t}) - Sl{t}) - Sl{t}) 
+ [x2(t)Ax2, b2)(pn{x3, 62) + (2 + a;2 + X3)0^(x2, 62)0^3^3, 

X as{t})hf\x,,X2,Xs,b,,b2)]exp[-SBit}) - SUt}) - SUt})]} . (17) 
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Note that when doing the above integrations over Xi and 6,, we have to include the correspond- 
ing Wilson coefficients Ci evaluated at the corresponding scale tj. The expression of Wilson 
coefficients are channel dependent which are shown later in this section. The functions hi, 
coming from the Fourier transform of H, are given in Appendix In the above equations, we 
have used the assumption that xi -C X2,xs. Since the light quark momentum fraction xi in 
B meson is peaked at the small region, while quark momentum fraction X2 of pion is peaked 
at 0.5, this is not a bad approximation. After using this approximation, all the diagrams are 
functions of fcf = Xiuib/V^ of B meson only, independent of the variable of kf. For example, 
by calculating the diagrams (b) we shall demonstrate it. 

(7r(P2)7r(P3)|0^^|S(pB)) 

d%d% Mhmh) ^2(p. _ kt)k:,%i]l2{kt - kt)k^ + 

d%d% { (,p.,,^,tM^,.) + ^ (^) ) ' ^''^ 

where the momenta are assigned in Figure |^. The calculation from second formula to last one 
is approximated as (ki) -C {k2)- This approximation is equal to take the momenta of spectator 
quark in the B meson as ki = (0, fcf , ki^). We neglect the last term which is higher order one 
in terms of I/itlb expansion. Therefore the integration of eq.@ is performed safely. Though we 
calculated the above factorization formulae by one order in terms of ag, the radiative corrections 
at the next order would emerge in forms of a^ln(m/t), where m's denote some scales, i.e. m^, 
1/6, . . ., in the hard part H{t). Selecting t as the largest scale in m's, the largest logarithm 
in the next order corrections is killed. Accordingly, the scale tj's in the above equations are 
chosen as 

tl = max{y/x^mB,l/bi,l/b2) , (19) 
tg = max{y^mB, l/6i, I/&2) , (20) 
td = meLx{y/xiX2mB, y/x^niB, l/h, I/62) , 

= max(./r^mB, 1/61, 1/62) , 
tj = meix{y/x^mB, \/x2 + - X2X3mB, l/h, I/62) , 
ta = max{y/x^mB, I/62, l/h) ■ (21) 



They are given the maximum values of the scales appeared in each diagram. 

In the language of the above matrix elements for different diagrams eg. ( P"3| - |T7D , the decay 
amplitude for tt+tt^ can be written as 



fnFe 



a kCi + C2 - 6 + -C3 + Cio + -Cc 



+ Me[^uCi-UC3 + Cg)] 
+ Ma 

— fsFait 



iuC2 — ^t[C^ + 2C4 + 2Cq + -Cs — -Cg + 2^10 



-C5 + — -C7 — -Cs 
6 b 2 



where ^„ = V^*^^^; 6 = ^tb^td- The decay width is expressed as 



1287r ' ' 



(22) 



(23) 



The C'^s should be calculated at the appropriate scale ti using eg. ( |6T| , |62D in the appendices. 
The decay amplitude of the charge conjugate decay channel n^n~ is the same as eg.(|22|) 

except replacing the CKM matrix elements to ^* and to Q under the phase convention 
CP\B°) = 15°). 



The decay amplitude for B^ 
-V2MiB^^n\^) = 



0^0 



TT TT 



can be written as 



fnFe 



'1 3 15 

+6 ( gC's + C4 + -C7 + -Cs — -Cg — Cio 



+ M, 

- Ma 

+ fBFa^t 



^uC2 — ^Cs + 2C4 + 2C6 + -Cs — -Cg + 2^10 
-C5 + Cg — -C7 — -Cg 

3 D / 



(24) 



The decay amplitude for B^ tt+tt^ can be written as 

■4 



a(Ci + C2) - 6 (2C10 + 2C9 - ^C7 - ^Cs 

\C, + ^Cy 



+ M, 



e«(Ci + C2)--6(C8 + C9 + Cio) 



(25) 
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From the above equations (|22| , p^ , |25| ) , it is easy to see that we have the exact Isospin relation 
for the three decays: 



(26) 



4 Numerical calculations and discussions of Results 



In the numerical calculations we use 



Agr^^ = 0.25 GeV, = 0.13 GeV, Jb = 0.19 GeV, 



Mb = 5.2792 GeV, Mw = 80.41 GeV, 



tb± = 1.65 X 10-^^ s, tbo = 1.56 x 10 



-12 



(27) 



and 



4.5 MeV, rud = 1.8m„, (28) 

which is relevant to taking mo = 1.5 GeV. For the vr wave function, we neglect the b dependence 
part, which is not important in numerical analysis. We use 

3 



(f)^{x) 



:/^x(l -x)[l + a^(5(l - 2xy - 1)], 



(29) 



with = 0.8, which is close to the Chernyak-Zhitnitsky (CZ) wave function |2^. For this 
axial vector wave function the asymptotic wave function a"^ ~ , is suggested from QCD 



sum rules , diffractive dissociation of high momentum pions |^ , the instanton model ||24 
and pion distribution functions [E^], etc., but we adopt = 0.8 according to the discussion 



in ref. |^6|. 0^ is chosen as asymptotic wave function 

3 



Ux{l-x)[l + a%5{l-2xy-l)], 



(30) 



with = 0. For B meson, the wave function is chosen as 
(f)B{x,b) = Nbx'^{1 — xY exp 



Ml x^ 1, ,,2' 



(31) 



with uJbi = ujfy2 = 0.4 GeV [^, and A'^^ = 91.745 GeV is the normalization constant. In 
this work, we set uJbi = 001,2 for simplicity. We would like to point out that the choice of the 
meson wave functions as in eqs. (^^) and the above parameters can not only explain the 
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experimental data oi B ^ nn, but also B Kn p8| , pq| , Dn etc., which is the result of a 
global fitting. However, since the predicted branching ratio of B ^ mi is sensitive to the input 
parameters fs, mQ, a"^, and lj^i, we will at first give the numerical results with the above 
parameters, then we give the allowed parameter regions of fs, ttiq, a^, and Ubi constrained 
by the experimental data of i? — > vr^vr^ presented by CLEO. 

The diagrams (a) and (b) in Fig.^, calculated in eq.(|l3|) correspond to the S — > tt transition 
form factor F^'^{q'^ = 0), where q = Pb — -^2- Our result is F^'^{0) = 0.25 to be consistent with 
QCD sum rule one. This implies that PQCD can explain the transition form factor in the B 
meson decays, which is different with the conclusion in ref . . In that paper, because mo was 
not considered, perturbative contributions to F^'^(O) were predicted to be much smaller than 
nonperturbative ones. 

Although we take the CZ like wave function {a^ = 0.8) for 0^, one finds that the above 
parameters give the pion electromagnetic form factor to be consistent with the experimental 
data. The pion electromagnetic form factor in PQCD is given as 

FiriQ^) = IQttCf dx2dx3 b2db2b3db3 as{t) /ie(a;3, 0:2, 63, ^2) 
JO Jo 

x[x2Q'^(p7T{^2,b2)4>n{x3,b3) + 2ml{l - X2)4>'^{x2,b2)4>'^{x3,b3)^ 

X eM-Slit) - Slit)], (32) 

where —Q^ is the momentum transfer in this system, the scale t is chosen as t = max(^/a?2Q, I/&2, V^s) 
and m^'s are replaced by Q in the he, Si and 5*^. One may suspect that around xi,X2 ~ 0, 
the gluon and virtual quark propagators give rise to IR divergences which can not be canceled 
by the wave functions. However, in PQCD, the transverse momenta k^p save perturbative cal- 
culations from the singularities around a;i 2 ~ 0. There are still IR divergences around ~ 0, 
but the Sudakov factor which can be calculated from QCD corrections does suppress such a 
region, i.e., non-perturbative contributions, sufficiently. We show the dependence of Ft^^Q"^) 
(eq.(P^)) in Figure ^ with the experimental data |^. This figure shows that the parameters 
we used don't conflict with the data. We also show F.„{Q^) for = 0.8, 0.4, and 0. It indicates 
that Fj^iQ"^) is fairly insensitive to a^. 
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The CKM parameters we used here are 

\Vud\ = 0.9740 ±0.0010, 
\VJ = 0.0395 ±0.0017, 



IWKbl = 0.08±0.02, 
\VtlVtd\ = 0.0084 ±0.0018. 



We leave the CKM angle 02 as a free parameter. 02's definition is 
In this parameterization, the decay amplitude of S ^ vrvr can be written as 



(33) 



(34) 



M 



v:,VudT - v:,vuP 



tb vtdi 



1 + ze 



i{<t>2+S) 



(35) 



where z 



V^Vud 



= V:,VudT 

, and 5 is the relative strong phase between tree (T) diagrams and penguin 
diagrams (P). z and 5 can be calculated from PQCD. For example, in tt^ti' decay, we 

get z = 30%, and 6 = 130°, if we use the above parameters. Here in PQCD approach, the 
strong phases come from the non-factorizable diagrams and annihilation type diagrams (see (c) 
~ (h) in Figure |^). The internal quarks and gluons can be on mass shell providing the strong 
phases. This can also be seen from eq. (|^^-|BD|) , where the the modified Bessel function KQ^—if) 
has imaginary part. Numerical analysis also shows that the main contribution to the relative 
strong phase 5 comes from the annihilation diagrams, (g) and (h) in Figure |^. From the figure, 
we can see that they are factorizable diagrams. B meson annihilates to qq quark pair then 
decays to tttt final states. The intermediate qq quark pair represent a number of resonance 
states, which implies final state interaction. In perturbative calculations, the two quark lines 
can be cut providing the imaginary part. The importance of these diagrams also makes the 
contribution of penguin diagrams more important than previously expected. 

This mechanism of producing CP violation strong phase is very different from the so-called 
Bander-Silverman-Soni (BSS) mechanism |^, where the strong phase comes from the per- 
turbative penguin diagrams. The contribution of BSS mechanism to the direct CP violation 
in i? ^ TT+TT^ is only in the order of few percent [^0. It is higher order corrections (a^ 
suppressed) in our PQCD approach. Therefore in our approach we can safely neglect this 
contribution. The corresponding charge conjugate B decay is 



M 



VubV:^T - Vt,v;,p 



VubV:,T 



tbVtd-^ 



1 ± ze 



i{-<i>2+sy 



(36) 
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Therefore the averaged branching ratio for B ^ tttt is 

Br = {\M\^ + \M\^)/2 



= KbV^aTf [1 + 2z cos 02 cos 6 + z'\. (37) 

From this equation, we know that the averaged branching ratio is a function of CKM angle 02, 
if zcos6 7^ 0. 

The averaged branching ratio of B^ —>■ tt+tt^ decay which is predicted from the formulae in 
the previous section is shown as a function of 02 in Figure |^. To consider tuq required from chiral 
symmetry is essentially different with previous paper [Tl[|. This figure shows that mg enhances 



the branching ratio to agree with the experimental data. There is a significant dependence 
on the CKM angle 02. The branching ratio of B^ tt+tt^ is larger when 02 is larger. The 
reason is that the penguin contribution is not small. The CLEO measured branching ratio of 

B TT + TT- 

Br{B -> TT+TT-) = {4.3t\i ± 0.5) x 10-^ (38) 

is in good agreement with our predictions. This prefer a lower value of 02. However, the 
predicted branching ratio is sensitive to the parameters of input. Especially it is sensitive to 
Jb, rnQ and the meson wave functions. Therefore, it is unlikely to use this single channel to 
determine the CKM angle 02. 

The branching ratios of -B — tttt are sensitive to some input parameters. We give the 



parameter regions allowed by the experimental data in eq. (|38D. Rerevant parameter are mo, 
, and Uhi. Others are specified in the begining of this section. Here we check the sensitivity 
of our calculation on parameter mo, , and Uhi- First we fix mo = 1.5 GeV and show the 
allowed region for and uJbi- This is shown in Figure ^(a). One finds that the branching ratio 
is fairly insensitive to . Second we fix = and show the allowed region for Ubi and mo. 
This is shown in Figure ^(b). We see that the allowed region for and mo is quite large. The 
dependence on a"^ for the branching ratio of 5 — vr+vr" is given in Figure ^ As discussed in 
ref. the central value of the experimental data Rd = Bi{B^ D^7t~)/'By(B^^ D^n^) 



requires = 0.8, but this figure indicates that B n^n^ decay mode gives no significant 
restriction on a-^. Therefore, these figures show that the above set of parameters we choose for 
Figure ^ is in the allowed region, and that parameter space producing the experimental data, 
eq. (P^D, is quite large. 
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The branching ratio of —>■ vr+vr'^ has httle dependence on (1)2. It is easy to understand 
since there is only one dominant contribution from tree diagrams. The QCD penguin contribu- 
tion is canceled by isospin relation and the electroweak contribution is very small giving only a 
slight dependence on 02- The branching ratio of this decay is predicted as 3 x 10~®, using the 
parameters we list in the beginning of this section. 

For the decay of B the situation is similar to that of B^ ti^ti . There are 

large contributions from both tree and penguin diagrams. We show the averaged branching 
ratio of B^ -k^-k^ as a function of 02 in Figure ||. Although the branching ratio is small, the 
dependence of 02 is significant. The predicted branching ratio of B^ -k^-k^ is less than 10^^. 
This is difficult for the B factories to measure the separate branching ratios of B^ and B^ . In 



this case, the proposed isospin method to measure the CKM angle 02 P3 does not work in the 



B factories, since it requires the measurement of B^ tt^tt^ and B^ tt^tt^. 
Using eg. (|35|J36D , the direct CP violating parameter is 

.... _ I-MP-W 



^CP 



\M\^+\M\^ 

—2z sin 02 sin 6 



(39) 



1 + 2z cos 02 COS S + z'^ 

The direct CP asymmetry is nearly proportional to sin 02. We show the direct CP violation 
parameters (percentage) as a function of 02 in figure |^. Unlike the averaged branching ratios, 
the predicted CP violation in B decays does not depend much on the wave functions. They 
cancel each between the charge conjugate states shown in the above equation. The direct CP 
violation parameter of B^ tt+tt^ and tt^tt^ can be very large, which can be as large as 
40%, and 20% when 02 is near 70°. Because there is no annihilation diagram contribution 
in B^ Tr+Tr*^, the penguin contribution is negligible. The direct CP violation parameter of 
B^ — s> Tr+Tr*^ is also very small. It is a horizontal line in Figure ^ 

For the neutral 5° decays, there is more complication from the B^-B^ mixing. The CP 
asymmetry is time dependent [|, [3^ : 



Acpit) ~ cos(Amt) + a,+,f sin(Amt), (40) 

where Am is the mass difference of the two mass eigenstates of neutral B mesons. The direct CP 
violation parameter Af^p is already defined in eg. (|39D . While the mixing- related CP violation 
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parameter is defined as 

-2Im{Xcp) 
- 1 + |AcpP ^ 

where 



(41) 



V,lVu{f\H^ff\B') 
Using equations (^,^), we can derive as 

Usually, people believe that the penguin diagram contribution is suppressed comparing with the 
tree contribution, i.e. 2; -C 1. Such that Xcp — exp[2i02], cte+e' = — sin 202, and A^p ^ 0. That 
is the previous idea of extracting sin 202 from the CP measurement of tt+tt". However, 

z is not very small. From Figure |TD|, we can see that a^+e' is not a simple — sin 202 behavior 
due to the so called penguin pollution. 

If we integrate the time variable t, we will get the total CP asymmetry as 

ACP = T^^^P + T^«e+.', (44) 
1 + x"^ 1 + x"^ 

with X = Am/F ~ 0.723 for the B°-B° mixing in SM [H. 

The integrated CP asymmetries of B^ tt+tt^ and B^ tt^tt'' are shown in Figure |ll. 

Unlike the averaged branching ratios, the CP asymmetry is not sensitive to the wave functions. 



since these parameter dependences canceled out. It is rather stable. If we can measure the 
integrated CP asymmetry from the experiments, then we can use this figure to determine the 
value of 02. 



5 Summary 

We performed the calculations of B'^ tt^tt^, B^ — > tt^tt^, and B^ —* decays, in a 

perturbative QCD approach. In this approach, we calculate the non-factorizable contribu- 
tions and annihilation type contributions in addition to the usual factorizable contributions. 
The predicted branching ratios of B^ —>■ tt^tt^ is in good agreement with the experimental 
measurement by the CLEO Collaboration. 

We found that the annihilation contributions were not as small as expected in a simple 
argument. The annihilation diagram, which provides the dominant strong phases, plays an 
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important role in the CP violation asymmetries. We expect large direct CP asymmetries in 
the decay of tt+tt^, and ti^ti^. The ordinary method of measuring the CKM angle 

02 will suffer from the large penguin pollution. The isospin method does not help, since the B 
factories can not measure well the small branching ratio of B^ tt^tt^. Working in our PQCD 
approach, we give the predicted dependence of CP asymmetry on CKM angle 02. Using this 
dependence, the current running B factories in KEK and SLAC will be able to measure the 
CKM angle 02- 
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A Wilson Coefficients 



In this appendix we present the weak effective Hamiltonian Tiejf which we used to calculate the 
hard part H{t) in eq.(0). The Tieff for the AB = 1 transitions at the scale smaller than mw is 
given as 



n 



Gp 

71 



10 



We specify below the operators in TYg/f for b ^ d: 



(45) 



Oi = da'y^LUfS ■ Ufi-f^Lba , O2 = daYLua ■ Up'^i^Lhp , 

O7 = \dal^Lha ■ J2q' eq'q'p'Jf.Rq'fs , Og = ^d^^YLbp ■ Eg' eq^q'^-ff^Rq'^ , 
O9 = ^da'J^Lba ■ Eg' dq'q'plixLq'p , 0^ = ^da^Lbp ■ Eg' (iq'q'ulpiLq'^ ■ 

Here a and (3 are the SU{'i) color indices; L and R are the left- and right-handed projection 
operators with L = (1 — 75), R = {1 + 75). The sum over q' runs over the quark fields that are 
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active at the scale ^ = 0{mb), i.e., {q'e{u, d, s, c, b}). 

The PQCD approach works well for the leading twist approximation and leading double 
logarithm summation. For the Wilson coefficients, we will also use the leading logarithm 
summation for the QCD corrections, although the next-to-leading order calculations already 
exist in the literature ||T^. This is the consistent way to cancel the explicit /i dependence in 
the theoretical formulae. 

At mw scale, the Wilson coefficients are evaluated for leading order as: 

C2W = 1; 

a,„ = 0, i = 1, 8, 10, 



247r 671 sin^ 9 



w 



^Aw — 

Stt 

Stt 

= ^(4Co + /^o) 

OTT 



9w 



a 



ACo + Do + -4^(105o - 4Co) 
sm Uw 



where 



(47) 



^ I ( X a; , \ 

= 4ir^ + (^3T)^^""j' 

4 , -19x3 + 25x2 x2(5x2 - 2x - 6) , 

Uq = mxH ; 1 ; -; lu X, 

" 9 36(x-l)3 18(x-l)4 

2, x(x2 + llx-18) x2(4x2 - 16x+ 15) , 

Eo = — lnx+^ — - + — z ^Inx, (48) 

" 3 12(x-l)3 6(x-l)4 ' ^ ^ 

with X = m1/m^r. 

If the scale mi, < t < mw, then we evaluate the Wilson coefficients at t scale using leading 
logarithm running equations (^), in the appendix |C|. In numerical calculations, we use as = 
Ait/[(3i hi{t'^ / AqIj^ )] which is the leading order expression with A.qI<^ = 193MeV, derived from 
^QCD ~ 250MeV. Here jSi = (33 — 2nj)/3, with the appropriate number of active quarks Uf. 
Uf = 5 when scale t is larger than mj,. 
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The Wilson coefficients evaluated at t = = 4.8GeV scale using the above equations are 





= -0.27034, 


C2 = 


1.11879, 


C-s 


= 0.01261, 


C,= 


-0.02695, 


Cs 


= 0.00847, 


Ce = 


-0.03260, 




= 0.00109, 


Cs = 


0.00040, 


c. 


= -0.00895, 


Cio = 


= 0.00216. 



If the scale t < mi,, then we evaluate the Wilson coefficients at t scale using the input 
of eq.(^), and the formulae in appendix for four active quarks (rij = 4) (again in leading 
logarithm approximation) . 



B Formulas for the hard part calculations 

In this appendix we present the explicit expression of the formulas we used in section 3. First, 
we show the exponent s{k, b) appearing in eq. ( p3H55| ). It is given, in terms of the variables. 



q = \n{k/A), b = \n{l/bA) 



(50) 



by 



s{k,b) 



2 



gin 



A(3l 



q + b 



-In 21 



- 1 



In f 



(51) 



The above coefficients Pi and A^"^^ are 

33 - 2nf 



Pi 



12 
67 



,,2) 67 10 8^ , fe'^\ 

A*^^' = Uf + -/?iln 

9 3 27 ^ 3^ V 2 / 



(52) 



where 'Je is the Euler constant. 

Note that s is defined for q > b, and set to zero for q < b. As a similar treatment, the 
complete Sudakov factor exp(— 5) is set to unity, if exp(— 5) > 1, in the numerical analysis. 
This corresponds to a truncation at large /cy, which spoils the on-shell requirement for the light 
valence quarks. The quark lines with large kj- should be absorbed into the hard scattering 
amplitude, instead of the wave functions. 
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^-SB(t) ^ ^-sl{t) ^ Q~sl{t) ^gg^ ^]-^g amplitudes are expressions abbreviated to combine 
the Sudakov factor and single ultraviolet logarithms associated with the B and vr meson wave 
functions. The exponents are defined as 

Ssit) = 4^^^B/V2,h)-j^\n^^l^y (53) 
Slit) = s[x2mB/V2,b2)+s[{l-X2)mB/V2,b2)-j^\n^^l^, (54) 

Slit) = s {x,mn/V2, 63) + ^ ((1 - Xs)mn/V2, h) - 1 In • (55) 

The last term of each equations is the integration result of the last term in eq.(0). 
The function /ij's, coming from the Fourier transform of hard part H, are given as, 

he{xi,X2MM) = Kq {^JxiX2mBhi) [d{hi - b2)Ko (v^"^b&i) (v^"^b^2) 

+6{b2 - bi)Ko {^/x>2mBh2) h {y/x^niBhi)] , (56) 
hd{xi, X2, X3, bi, 62) = Kq i-iy/x2XsmBb2) [6'(6i - 62)^0 ( Va^mB&i) Iq (y/x^niBbi) 

+9{b2 - bi)Ko {^JxiX2mBb2) Iq {^/x^mBbl)] , (57) 
hf\xi,X2,X3,bi,b2) = Koi-iy/x^niBbi) 

X [9{bi - b2)Ko {-i^/x^niBbi) Jo (^2^^562) 
+6{b2 - bi)KQ {-iy/x^mBb2) Jq (y/x^niBbi)] , (58) 
hf\xi,X2,X3,bi,b2) = Ko (yxT+a^T^^x^msfei) 

X [9{bi - b2)Ko {-iy/x^TJlBbi) Jq {y/x^mBb2) 

+9{b2 - bi)Ko {-i^/x^mBb2) Jq {^/x^rriBbi)] , (59) 
ha{x2, X3, 62, ^3) = Kq {-i^/x^niBb^) [d{b2 - b3)Ko {-i^/x^mBb2) Jq {^/x^mBbs) 

+ 0{b3 - 62)^0 (-iy/x^mBbs) Jq (v^"^b&2)] , (60) 

with Jo the Bessel function and Kq, Iq modified Bessel functions Kq^—ix) = —{'k/2)Yo{x) + 
tin/2)Jo{x). 



C Wilson coefficients running equations above scale 

In this appendix, we list formulas for renormalization group running from rriw scale to t scale, 
where t > nib. These formulas are derived from the leading logarithm QCD corrections with 
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.1456 



0.8994 



active quarks |jl6 

Ci = i _ ^12/23 

Cs = 0.05107/-°-^°^'^ - 0.0714r/"'^/23 + 0.0054r/-°-^^^^ 
-0.1403r/0-^23° - 0.0113r/°-«99^ + l/eT]^^/^^ 

+C3^(0.286877-°-^°«6 + 0.0491r/-o i^^6 + 0.657^-^230 ^ o.OOGlr^o-ss^^) 
+C4^(0.3287r/-°-^°«6 + 0.0424r7-°-^^^'^ - 0.3263r7°-^23o _ o.0448r7°-«99^) 
+C5^(-0.0629r/-™ ^ o.l629r/-°-i456 - 0.1846r/0-^230 ^ o.0846r7°-«99^ 
+a»(0.0447r7'°-^°«6 - 0mQ3r]~^-'^'^ - 0.2610r7°-^230 ^ o.2226r/°-«99^) 

+C9^ (-0.032577-°-^°^^ + 0.0357r/"^/2^ - 0.0016r/-° 
+ 0.2342770-4230 _ o.25r/i2/23 ^ o.0141r/0-8994 

+C7^ (-0.00637]-°-^°^^ + 0.0163//-°-^'^^ - 0.0185r/°-'230 ^ 0.00857^' 
C4 = 0.0984//-°-^°^'^ - 0.071477''^/23 ^ o.0026r/^°-^^^^ 
+0.1214770-^230 _ 1/6^12/23 ^ o.OlSer^o-sss^ 

+C3^(0.5539r7-o-4086 ^ o.0239r/-0-i^56 _ o.5693r/0-4230 _ o.0085r/0-«994) 
+C4^(0.63 4877-°-^o«6 + 0.0206r7-°-i^5'^ + 0.2823r/°-^230 ^ o.0623r/°-«99^) 
+C5^(-0.1215r/-™ + 0.0793r7-°-i45^ + 0.159777°-4230 _ o.nrSr/^-^^^^ 
+C6^(0.0864r/-°-^°«6 - 0.0031r/-°-i^56 + 0.2259770-^230 _ o.3092r7°-«99^) 
+C9^ (-0.06277]-°-^°^^ + 0.0357r/~*^/23 _ 0.0008//"° 

- 0.2027r/°-^230 ^ 0.2577^2/23 _ o.0196r/°-«99^ 
+C7^ (-0.012277-°-^°^^ + 0.00 797/-°-^'^^ + 0.0160r7°-'230 _ 0.011777' 
C, = -0.03977/-™ + 0.0304r7~°-i^^6 + 0.01 17r7°-^230 _QQQ25^o.8994 

+C3^(-0.2233r7-o-^°«6 + 0.2767r]-'-"''^^ - 0.0547r70-^230 ^ o.0013r/°-«994 
+C4^(-0.2559r/-°-^°«^ + 0.2385r/-°-i^56 ^ o.0271r7°-4230 _ o.0098r7°-«99^ 
+^5^(0.049077-°-^°^^ + 0.917177-°-!^^'^ + 0.01 5477^-^230 ^ o.018577°-899^) 
+C6^(-0.03 4877-°-^°«^ - 0.035777-°-!^^^ + 0.0217r/°-^230 _^ o.048877°-«99^ 



■.1456 



0.8994 
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c« = 



+C9^ (0.0253r/-°-^0«6 - 0.0089r/~o-i^56 - 0.0195r/°-^230 ^ o.0031r7°-«994 
+Cj^ (0.0049r7-°-^°«6 

+ 0.0917r7-°-i^56 - 0.1r7-3/23 ^ o.0015r7°-^230 ^ o.OOWr/"-^^^^) , 
0.0335r/-°-^o«6 - 0.0112r/-° i^56 + 0.0239r7°-^2=^° - 0.0462r70-8994 
+C3^(0.1885r7-o-^o«6 - O.lOlTr^-^-i^^^ - 0.1120r7°-^230 ^ o.0251r/°-«99^) 
+C4^(0.2160r/-°-^o«6 - 0.0877r/-°-i^56 + 0.0555r7°-4230 _ o.i839r7°-«99^) 
+C5^(-0.0414r/-°-^°«6 - 0.33 7077-°-i^^^ + 0.0314r7°-^230 ^ o.3469r/°-«99^ 
+a»(0.0294r7-°-4°8'5 + 0.0131r/-0-i^56 ^ 0.0444770-^230 ^ o.9131r7°-«99^) 
+C9^ (-0.021377-°-^°«^ + 0.0033r/-°-i45^ - 0.0399r/°-^230 ^ o.0579r7°-«99^ 
+^7^ (-0.0041r/-°-^°«6 - 0.0337r/-°-i^56 + 77-2/2730 + 0.0031r70-^230 
+ 0.0347770-899^ - 7/24/23 /30) , 
C7 = C7^r7-='/'^ 
C, = ^C7. (-77-2/^2 + 7724/^2) , 



= iCg^ (77-6/23 +77 

1 



12/23 



ClO = iCg. (77-6/23 - ,712/23) ^ (61) 

■where 77 = as(t)/as{mw)- 

D Wilson coefficients running equations below scale 

In this appendix, we hst formulas for renormahzation group running from mj, scale to t scale, 
-where t < rrif,. These formulas are derived from the leading logarithm QCD corrections with 
four active quarks |]16 . 

CC^ = C4 (0.3606C-°-2^69 + o.03166C-°-^2^^ - 0.3626C°-^^°^ - 0.0297C°-8^^^) 

+C10 (0.0149C-°-=^''' - 0.0020C-°-'=^'' - 0.498lC°' + 0.5^'/'" - 0.0148^'') 
+C2 (0.0651C-°-2^69 _ 0.0833C-6/25 
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+ 0.0046C-°-'''' - 0.2265C°-''°' + 0.25C''/'' - 0.0099C°-'''') 
+C3 (0.3308C-°-'''' + 0.0356C-°-'''' + 0.6337C°-"^°^ - O.OOOIC^'""'') 
+Ci (0.0502C"°-^^^^ - 0.0833C"^/^^ + 0.0066C"°-^^^^ + 0.2717C°-^2°^ 
_Q 25^12/25 ^ o.0049C°-'''') + Cg (-0.0149C-°-'''' 
+0.0020C-°-'''' + 0.4981C°-"'°' - + 0.0148C°-'''') 

+C5 (-0.0598C-°-'''' + 0.1371C-°-'''' - 0.1473C°-''°' + 0.0700C°-'''') 
+C6 (0.0377C-°-'''' - 0.0045C-°-'''' - 0.2210C'^-''°' + 0.18775C°-'^'') 
+C7 (-0.0150C-°-=^''' + 0.0343C-°-'''' - 0.0368C°-'20i ^ o.0175C°-'''') 
+Cs (0.009C-°-'''' - 0.0011C-°-'''' - 0.0553C°-''°' + 0.0469C°-'''') , 
CC4 = Ce (0.0640C-°-'''' - 0.0021C-°-'''' + 0.2018C°-''°' - 0.2637C°-«^5^) 

+C5 (-0.10156C-°-''*^' + 0.06538C-°-'''' + 0.1345C°-''°' - 0.09836C°-'''') 

+C9 (-0.02528C"°-^^^^ + 0.0009C"°-^^^^ 

- 0.4549^°' + - 0.0207C'') 

+Ci (0.08515C"°-^^'^^ - 0.0833C"^/^^ + 0.0031C"°-^^^^ 

- 0.24809C° ''°' + 0.25C''/'' - 0.00688C° '^'') 
+C3 (0.5615C-°-'''' + 0.01699C-°-'''' - 0.5787C°-"^'^^ + 0.0002C° «'5i) 
+C2 (0.1104C"°-^^^^ - 0.0833C"^/^^ 

+0.0022C-'^-'''' + 0.2068C°-^2°^ - 0.25C''/'' + 0.0139C"-'^'') 

+C10 (0.0253C"°-^^^^ - 0.0009C"°-^^^^ 

+ 0.4549C°-"'°' - + 0.0207C°-'''^) 

+C4 (0.6121C-°-'''' + O.OISIC""-'"'" + 0.3311C°-"^°^ + 0.0417C°-'"'') 
+Cs (0.0160C-°-=^''' - 0.0005C-°-^'^' + 0.0505C°-''°^ - 0.0659C°-'''^) 
+C7 (-0.0254C-°-'^''^ + 0.0163C-°-'''' + 0.0336C°-''°' - 0.0246C°-'''') , 
CC5 = C4 (-0.2291C-°-'"'' + 0.2167C-°-'''' + 0.0192C°-''°' - 0.0067C°-«^5^) 

+C10 (-0.0095C-°-'''' - 0.0136C-°-'''' + 0.0264C°-''°^ - 0.0034C°-'''^) 
+C2 (-0.0413C-'^-'^''^ + 0.0316C-°-'''' + 0.0120C°-''°' - 0.0022C° «'5i) 
+Cs (-0.2102C-°-''^' + 0.2438C-°-'''' - 0.0336C°-''°') 
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+Ci (-0.0319C-°-'''' + 0.0452C-°-^'^' - 0.0144C°-''°^ + 0.0011C°-'''') 
+C9 (0.0095C-°-'''' + 0.0136C-°-'''' - 0.0264C°-''°' + 0.0034C°-'"'') 
+C5 (0.0380C-°-=^''' + 0.9382C-°-^-'^' + 0.0078C°-''°' + 0.0159C°-'''') 
+C6 (-0.0240C-°-'''' - 0.0305C-°-'''' + 0.0117C°-''°' + 0.0427C°-'''^) 
+Cs (-O.OOeOC-^-"""'^ - 0.0076C-°-'''' + 0.0029C°-''°' + 0.0107C°-'''^) 
+C7 (0.0095C"°-^^^^ + 0.2346C"°-^^^^ 
- 0.25C-=^/'' + 0.0020C°-^2'^^ + 0.0040C°-'^'') , 
CCq = C4 (0.1825C~°-^^^^ - 0.0784C~°-^^^^ + 0.0449C°-'^^°^ - 0.14894C°-^^^^) 

+C10 (0.0075C-°-'''' + 0.0049C-°-'''' + 0.0617C°-''°^ - 0.07412C°-'''^) 
+C2 (0.0329C-°-'''' - 0.0114C-°-'''' + 0.0280C°-''°' - 0.0495C°-'''') 
+C3 (0.1674C-°-'''' - 0.0882C-°-^'^' - 0.0784C°-''°^ - 0.0007C°-'''^) 
+Ci (0.0254C-°-'''' - 0.0163C-°-'=''' - 0.0336C°-''°' + 0.0246C°-'''') 
+C9 (-0.0075C-'^-'"''^ - 0.0049C-°-'''' - 0.0617C°-''°' + 0.07412C°-'''^) 
+C, (-0.0303C-°-='"'' - 0.3395C-°-'''' + 0.0182C°-''°^ + 0.35157C°-'''^) 
+Ce (0.0191C-°-'''' + O.OllOC-"-'"'" + 0.0274C°-^2°i + 0.94253C° «'^^) 
+Cs (0.0048C"°-^^^^ + 0.0028C~°-^^^^ + 0.0068C°-^^°^ 

+ 0.2356C°-'''' - 0.25C''/'') + (-0.0076C-°-'''' - 0.0849C-°-'''' 
+0.0833C~^/^^ + 0.0046C°-^^°^ + 0.0879C°-^^^^ - 0.0833C^^/^^) , 

CCs = + C''/'')/3. + C8C''/'^ 

CC,o = C9 (0.5C-*^/'' - 0.5C^'/'') + Cio (O.SC"'/'' + 0.5C^'/'') , 
where C = a,(t)/a,(ms). Here A^^^^ = 250MeV. 
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Figure Captions 



Figure 1: One of the decay process which contributes to B ^ tttt decay, b quark decays to 
produce a fast moving u quark. In general, this quark and d quark are not lined up to form 
a pion. A gluon exchange is necessary in order that these quarks are lined up to form a pion. 
The part enclosed by dotted line describes the six-quark effective operator. 

Figure 2: b-Q dependence of e"**. Nonperturbative region on 6 ~ 0(A~^) is suppressed by this 
exponent. Since the pion wave function has two Sudakov factor accompanied with two light 
quarks, this suppression becomes much stronger. 

Figure 3: Diagrams contributing to the B tttt decays. The diagram (b) corresponds to 
Figure 0. 
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Figure 4: Q"^ dependence for F^^^Q^) with the data The sohd, dashed, and dotted hnes 
correspond to = 0.8, 0.4, and 0, respectively. 

Figure 5: Averaged branching ratios (in unit of 10~^) of B^{B^) tt^tt", rriQ = 1.5 GeV 
(sohd hne), mo = GeV (dashed hne). The two dotted hnes indicate the la region of CLEO 
experiments in eq. 



Figure 6: Here we check the sensitivity of our calculation on parameter mo, , and uJbi- Others 



are defined in the begining of Section The shaded areas are allowed by the data, eq. (^ 
for arbitrary 02; a) we fix mo = 1.5 GeV and show the region for and Ubi] b) we fix = 
and show the allowed region for Ubi and mo. 

Figure 7: Dependence on a"^ for the branching ratio (in unit of 10~^) of B ^ tt+tt". = 
0.8, 0.4, in descending order, respectively 

Figure 8: Averaged branching ratios (in unit of 10^^) of B^ tt^tt*^ as a function of CKM 
angle 02- 

Figure 9: Direct CP violation parameters (in percentage) of B^ tt+tt^ (dotted line), and 
B^ tt+tt'^ (solid line), and B^ (dashed line) ClS db function of CKM angle 02- 

Figure 10: CP violation parameters a^+e' (in percentage) of B^ tt+tt^ (solid line), and 
B^ TT^n^ (dotted line), as a function of CKM angle 02- 

Figure 11: The integrated CP asymmetries (in percentage) of 5° tt+tt^ (solid line), and 
B^ — i> TT^TT^ (dotted line), as a function of CKM angle 02. 
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^"^'^575^ : hard gluon 

: quark line 

bound state to form a meson 



E°(at rest) 
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